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We study elasticity of a perpendicular flux lattice in a layered superconductor with Josephson 
coupling between layers. We flnd that for flux displacements p the energy contains In p terms, so 
that elastic constants cannot be strictly deflned. Instead we define effective elastic constants by a 
thermal average. The tilt moduli have terms ~ InT which for Aj <^ a, where Xj is the Josephson 
length and a is the flux line spacing, lead to (p^) ~ r/|lnT|. The expansion parameter indicates 
that the dominant low temperature phase transition is either layer decoupling at high fields (A,j ^ a) 
or melting at low fields (Aj <C a). 

bb"' 

0\ ■ I. INTRODUCTION 

a^ : 

The properties of Abrikosov flux lattices in layered superconductors are of considerable interest in view of numerous 
Ph ' experiments on high Tc compounds Q . For magnetic field perpendicular to the layers, the flux lattice can be considered 
as two dimensional point vortices in each superconducting layer which are stacked one on top of the other. Each point 
vortex, or a pancake vortex, represents a singularity of the superconducting order parameter, i.e. the superconductor's 
phase in a given layer changes by 27r around the vortex. These pancake vortices are coupled by their magnetic field as 
well as by the Josephson tunneling between nearest layers. The fluctuations of the displacements of pancake vortices 
I ' are manifested by a variety of experiments and affect phase transitions such as melting of the flux lattice and 
J>. , layer decoupling, i.e. vanishing of the interlayer Josephson coupling on long scales 

• A harmonic expansion of the flux lattice energy to second order in displacements of the flux lines, defines in 
0^ ', general the elastic constants For anisotropic layer systems a harmonic expansion was studied by Glazman and 
Koshelev (GK) This expansion is, however, nontrivial, since it involves expanding the nonlinear Josephson 
coupling J cos(j)„^n+i where (j)n,n+i is the relative phase of neighboring layers. When two pancake vortices which are 
one on the top of the other are separated by a distance 2p (parallel to the layers) then 4>n,n+i has large variations 
qI^ in a circle of radius p between these pancake vortices. This effect has led GK Q to anticipate a p^\a.p term in the 
"p-, , energy expansion when p > S, where ^ is the in-layer coherence length. This term was also found by Kramer ]lO| for 

■ a single vortex line and independently by us . 
^ In this work we present the detailed energy expansion. In addition to the difficulty at short scales, leading to 

I the p^ Inp term, we find that the convergence parameter of the expansion vanishes as 1/| ln_Ej| when the Josephson 
■ O coupling Ej 0. Thus the Ej = elastic constants are not recovered by a Ej limit. We then show (Section 
^ III D) how to define effective elastic constants and apply our results (Section IV) to thermal averages of (p^). In 
Q particular we find (p^) ~ T/|lnT| for Aj ^ a, where Xj is the Josephson length and a is the flux line spacing. 
• • ; The expansion parameter indicates the type of dominant fluctuation with a related instability, i.e. for Xj ^ a the 
. ^ ' decoupling transition dominates while for Xj ^ a the melting transition dominates. The relevance to experimental 
data is discussed in section V. 
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II. THE MODEL 



We start from the Lawrence-Doniach free energy in terms of superconducting phases (f>n{i^) on the n-th layer 
and the vector potential A{r, z): 



f d'rdz (V X Air, z)Y + A ^(^ V0„(r) - l(r, z)fS{ 



z — nd) 



d^rcos 0„(r) - 0„_i(r) - — / ^^(r,z)dz 

^0 J(n-l)d 



(1) 



where Xab is the penetration length parallel to the layers, d is the spacing between layers and $o = hc/2e is the flux 
quantum. 

In Appendix A we consider a simplified model with two layers in the e ^ (i.e. A{v,z) decouples from (/)ri(r)) 
limit and two pancake vortices, one at the bottom layer at Ri , and the other at the top layer at Ri + p. It is rather 
straightforward, as shown in the Appendix A, to show that the free energy expansion has a p^hip term. 
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We proceed with the partition sum for Eq. (||) which involves integrating over (/)„(r) and A(r, z), subject to a gauge 
condition. Since ^(r, z) is a Gaussian field (choosing the axial gauge Az{r, z) — 0, v4(r, z) = (A(r, z), 0) ) we can shift 
A ^ A + 5 A. where A(r, z) now satisfies the x, y components of 

V X V X A(r, ^) - Et?^-^" W - A(r, z)]S{z - nd) (2) 

ab n 

and then fluctuations in 5 A. decouple from those of 0n(r). The partition sum at temperature T is now 

Z = Jvcj)nir)eM-^/T] (3) 

with A{r,z) in Eq. given by the solution of Eq. (j^). Note that since Eq. (|^) is gauge invariant under A —>■ 
A — ■|^Vx(i", nd) and 0„(r) ^ 0n(r) — x(i": nd) one can in fact choose any gauge. 
We now decompose 0„(r) to 

,^„(r)-/„(r)+5^s„(r>(r-r') (4a) 

0„(r) = ^O(r)-0O_i(r) (4b) 

where (/>°(r) is the nonsingular part of 0n(r), a(r) = arctan(2//a:), s„(r) = 1 at pancake vortex sites and s„(r) — 
otherwise. The sum in Eq. (^) is then a sum on r' being the vortex positions on the n-th layer. 
Solving Eq. (|^) for A in terms of 6'„(r) and s„, substituting in Eq. (|I]) yields |^ 

T = T„ + Tj + Tf 

where T^, is the vortex- vortex interaction via the 3D magnetic field, J-j is interlayer Josephson coupling term and Tf 
is an energy due to fiuctuations of the nonsingular phase: 

= ^ s„(r)G'i,(r-r';n-n')s„'(r') (5a) 
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Tj^-EjY, j d2r(cos[r(r)+^(s„(r')-s„_i(r'))a(r-r')]-l) 



(5b) 



■^/^^EG7'('?'^)l^(q'^)l' (5c) 
Here (q, A;) is a 3D wave vector, z is the direction perpendicular to the plane of layers and 

d sinhgd 
4A^j,g smh^ ^ + sm^ ^ 



= l + -#sin^- (6c) 



For deviations u" of 2D vortices on the 7i-th layer from equilibrium positions R; of a hexagonal lattice , the function 
s„(r) is 

^ . ^ f 1 if r = + uj^ 
> otherwise. 

The Fourier transform 

u(q, k) — ^ u" exp(iqR; + iknd) 

n,l 

identifies longitudinal u'(q, k) = q - u(q, k)/q and transverse u*''(q, fc) = [q x z] • u(q, k)/q components of u(q, k) 
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u(q, k) = It' (q, k)q/q + u*'' (q, k) [q x z]/q. 

If the free energy can be expanded to second order in u(q, k) then the compression cn, shear cge and tilt C44 moduh 
for the vortex lattice are identified by: 

T=lJ^J^^'^{daY{[q'cMq,k) + klcUq,k)]\u\ci,kr + 

[q^Ceeiq,k) + klcUq,k)]\u'^ici,k)\^} (7) 

where k^ = sin^ a? is the area of a unit cell (a^ = '^q/ B) and /g^-i^ *1 integration over the Brillouin zone. 
We assume below that d <ti a, Xab as is the case for high compounds ^TJg] . 

Note that for q = there should be no distinction between transverse and longitudinal £44(0, k) = 0^4 (0, fc), however 
for q 0, 044(9, k) and 044(5, k) may differ. 

III. ELASTIC CONSTANTS OF THE FLUX LATTICE 
A. Magnetic Coupling 

We consider first the case with no Josephson coupling, Ej — 0. The vortex- vortex interaction has then the form 



The first and the second terms can be expanded with respect to q • u since in absence of the I — I' , n = n' term 
they converge. It is important not to decompose the (1 — Si^i') or (1 — Sn^n') factors until all integrals converge; the 
1 = 1', n = n' terms produce then the integral terms in the following elastic matrix 

= \ J^J^^^ ^^iday^''{n,k)uM,k)u*{q,k) 

<^'^"(q, fc) = (^)2 ^ [G„(|Q - q|, fc)(Q - q)»(Q - q)^' - a(g, 0)Q'g^"] 
" Q 

where Q are 2D reciprocal vectors of the hexagonal lattice. 

Considering q — > we use the symmetry of the hexagonal lattice 

i,j Q 

Y.siQ)Q^QJQl^o (8b) 

E 9{Q)QiQjQiQm = -^{SijSim + SiiSjm + SimSji)'^ g{Q)Q^ (8c) 



i.j,l.7n Q 



and separate the Q = and X^q^^^o parts. We consider flux line spacing a 3> fi, so that the Q sums involve many 
terms (~ (a/d)^) and the sums can be approximated by integrals, 

E « 752 / QdQ (9) 



Q#0 



Qo JQo 
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where Qq = 4ttB/^o = Att/o?, i.e. ttQq is the area of a Brillouin zone. 
Note that for Q : 



Gy{q, k) 



47r 



1 



+ 



1 



(10) 



of the tih moduh C44. Therefore the Ej — compression cn, shear cee and tilt C44, 044'^ moduU for q <^ Qo are 



The first term of Eq. ([10|) contributes to the compression moduh cn while the second term to the longitudinal part 
luh C44. 

cii(q, k) 



1 



2 'da 

3 

16^2 a 



QdQ^^ dQ' 



4 E 



1 B$o 



a^d^ 1 d 
(27r)2 p 9p 



[p^G,{p,k)] 



(27r)2 p dp dp 



[p-^—G,{p,k)] 



(11a) 



C66(g,fc) = (^)2^[^X]iA 

^ g$o 
(87rAafc)2 
52 1 



3 dGy{q,k) 

dQ ^ 



C44 {q,k) 



tr.O 



tr,Q, 



47r 1 + Xl,{q^ + fc2) 
1 



-M\l^k) - 1(^)2J- ^[a„(Q,fc) -a(Q,0)]Q2 



d^l d . ^dGy{p,k) 



(27r)2 p 9p 



[P 



dp 



2B$n 



Q#0 



[^^xi^^ki i+ek^. 



(lib) 

(Uc) 

(lid) 



Note that C44 (5, k) ^ C44 (g, fc) even for g ^ due to the singular form of the vortex- vortex interaction Gt,(q, fc) 
(Eq. (|o|)). At g = the C44 terms combine into ifc^ [044*^(0, fc) + 044(0, fc)]|u(0, fc)|2 which can also be verified by 
direct expansion for u(0, fc). As shown below, a finite Josephson coupling restores the equality 044(9, fc) = £44(5, fc) at 
a^O. 



B. Josephson Coupling: "Naive" Expansion 



We consider now the contribution of the Josephson coupling Eq. (5b) to the elastic constants by a conventional 
expansion, reproducing the results of GK |^ . 

The singular part of Josephson phase difference in the interlayer Josephson coupling term can be written as [see 
Fig.l] 

^r(r) = a{v - Rr - pD - a(r - + p^) (12) 

Here we defined 

Rr = Rz + + (13^) 

p7 = ' (13b) 

vr(r) = r-Rr. (13c) 
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FIG. 1. The Singular phase difference ip"ir) and the "p circle" where |r — R"| < p" . 



The usual way for treatment of the cosine term in Eq. ( pb| ) is by using a "naive" double expansion: 
i) Expansion of the cosine with respect to the phase difference 

ii) Expansion of the singular phase difference ipf'ir) with respect to pf 

^r(r)«-2Va(r-Rr)-pr 

Each pair of vortices, displaced by u" and u"~^ respectively, defines a "p-circle" in space where |v"(r)| < p", see 
Fig. 1. Within a "p-circle" ipri^) ti^s a 27r discontinuity and therefore cannot be expanded. The expansion (ii) is 
reasonable only in the region far from the "p-circle" where vf'{r) 3> |pf |. 

Within the approximations (i) and (ii) we can write 



E / d2r(cos[r(r)+E^r(r)]-l) 

n I 

^-^E / d'rr(r)+E'/'r(r)P 



I '^'^ -(|0(q, fc)|2 + r (q, A;)B(q, k) + c.c + |i?(q, k)\^) (14) 



2d J (27r) 



where we define 



S(q, /fc) = dE / d^re^^^^''"-' E ^^(r) = _l^!fciL^ 



i?r(r) = -2Va(r - RD • = ^^.^"/'^^^ (15) 



and use the Fourier transform 

2-!Ti[z X q] 



Note, that merely the use of the expansion (ii) leads to an error of order since the difference (V'['(r))^ — (_B"(r))^ 
is of order 1 in the "p-circle" with area ^ p^. 
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Combining Eq. (Q with of Eq. @ yields 



Here we introduced: 



El 

' 2d 



(fqdk 



_ -^iGjH^,k)+Ej/d)\e{q,k)-e^{q,k) \' + 
5 d'^qdk \4:nid[z x q] • p(q, fc)p 



(27r)3 



g2 +772 



(16) 



,7 V 2 



4A^ 



where the Josephson length is 



X.J - 



AXabVEjd^^ 



Since d <C Aafc typically rjk ~ S/Aj for most k averages below. 

The last term in Eq. contributes to the longitudinal C44 and transverse C44 part of the tilt moduli. Rewriting 
the integrand in the form 



/CXD I'Qo 
d^qg{q,k)=Y, d'q5(q + Q, fc) , 
Q 



(17) 



using the symmetry of the hexagonal latti ce Eqs. (^) as well as Eq. (^) with an upper cutoff l/i^, the tilt moduli 
[including the magnetic contribution, Eqs. (11c lid)] can be written as 



^44(91 k) — C44 



£44(9, k) — C44' + 



(8^Ac)2 
52 



1 2S$o 
4^TTW + Apf " (8^ 



lne^[Q2 + (i + A2,fc2)/^: 



(18a) 
(18b) 



where Ac = XabXj/d. The result for 0^4 (g, fc) was obtained by GK |]. 

It is interesting to note that for (7 — > and finite Josephson coupling Ej t he til t moduli are equal. Note also that 
the limits q ^ and Ej ^ do not commute in the second term of Eq. (18b). In fact we show below that the 
expansion of Josephson term breaks down when Ej 0. 

In summary, the "naive" expansion needs a revision to correct two aspects: (i) short scale behavior - expansion is 
not allowed in the "p circle", and (ii) long scale behavior - generation of an Ej independent term when q ^ 0. 



C. Josephson Coupling: Proper Expansion 



We proceed now to a method which avoids the " naive" expansion by expanding the cosine term in Eq. ( pq ) directly 
in terms of p. 

We define f?"(r) = ^"'-'^(r) + e"(r) and expand the Josephson coupling term in Eq. (5b) with respect to e"(r) to find 
an optimal (?"'^(r) for which the expansion is allowed: 



EjJ2 f dS cos[r'i(r)+^7Ar(r)]-l-ie"(r)'-e"(r)sinr^i(r)+^Vr(r)] 

n \_ I I 

+0(e^eV') 



(19) 
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We show below that terms of order e^(cos[0"'^(r) + ~ 1) contribute 0{e^p^) correction to the free energy 

after integration over r. 

The expansion is most efficient when the term hnear in e(q, k) vanishes. This determines 9^ (q, k) to be the solution 

of: 



(20) 



To solve Eq. (|20|) we introduce the functions 



^^(r) = e'[«r''('-)+'Ar(r)] I (21) 



Sr{r)^9-^\r) + Cr{r) (22) 
2[pY X vr(r)]. 



where C"(r) is defined as 



crir) 



w,"(r)2 + (p[')2 
Ciiq,k) =dJ2 I rf^re*''""+'''"''C;"(r) = 

n •' 

= AiTidY^ {qp\'Ki{qpY)) e^i^^r+^fc^d ^23) 

n ^ 

In the rest of this section we identify expansion parameters ( Eq. ( |26| ) below) which allow a solution of Eq. (|2^), and 
derive the free energy expansion (Eq. ( pl| ) below). 
Since the function C"(r) is close to 



2[pr X vrw 

[K"(r)' + |PrP)'-4(vr(r)-pr)' 



sin (r) = " "'^^^^^ r (24) 



for both <C and ^ p", the difference between imaginary part of D"(r), /mZ)"(r), and (5"(r) is only on the 
" p-circle" , so that 



d'^r\ImD]'{r) - (5,"(r)| - 0(p2) (25a) 

/ d'rJ2lmDi {r)ImDf,{r) = / d^r ^ 5['(r)5p(r) + 0(pV«)- (25b) 

1=^1' •' i=ii' 

We show now that an expansion in p" is possible if the following expansion parameters are small 
2d d^q r''^ 



Ori r^i'^ A^n ni<^ 
X=— ^ / dfc(|p*'^(q, k)\'} « 1 if Xj » a 



(26a) 

l/Aj </ J-tt/q 

(e> < 1 i/ Aj < a (26b) 

where (p^) is an average of which is diagonal in q, k. The case of thermal average is evaluated in section IV. The 
parameter x controls the expansion of the sine term in Eq. (20) and is evaluated in Appendix B while (e) — 0{p^) 
results from the solution of Eq. (|20| ) which is to leading order in p, so that the term linear in e in Eq. ( |l9|) survives 
and leads to higher order corrections. 

We claim then that the solution ^?"'^(r) = J^i ^z"'^!'^) ( Po| ) (compare with ^?"'"(r) from "naive expansion") is 



so that with Eq. (|22|) 

I d pco „3„m 



^r(r) = / * /°°,,'^vr^i(.Pryi(.<)[,„ , ^™]^,.(«-™). (28) 
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where Ki, Ji are conventional Bessel functions. The function i5"(r), in terms of v" = r — R", decays slowly as 
for a < w" < Xeff where Xj < Ae// < Ac depends on the configuration of p", but for u" > Ag// it decays 
as exp(— up/Ae//). The exponential decay allows the convergence of the I summations inside the sine in Eq. ([2(]); 
however, since the e xpon ential decay sets in at the scale Xeff which diverges when Ej — > 0, we expect the expansion 
parameter x of Eq. (26a) to diverge, i.e. the expansion is invalid when Ej — > 0. 

The convergence o Tv, SJ' jr) implies that sinQ^i SJ'ir)) can be expand ed. More precisely, as shown in Appendix 
B, the condition of Eq. (|26a|) leads to Eq. (B3), which together with Eq. (|25ci|) yields 



J d^rsin ^^^^^(r) + ^ Vz^W j 3''^'+''''"^ = 

n I •' 



Substituting in Eq. ( pO[ ) shows that Eq. ( ^7) is indeed the solution for Eq. 
Furthermore we have from Eq. (B4), using Eq. ( 25b| ) 



(29) 



I), i.e. it is the optimal 6 



n,l 



/dMcos(r^i(r)+^V^r(r))-l] = 
I 

Jd'r (^i?ei?r(r)-^E^'^^r(r)/mZ?P(r) I [1 
/ rf^r I ^ ReDUv) -\Y. Sr{r)6r,{r) | [1 + 0{x)] 



0{x)] = 



(30) 



where ReDf{r) is the real part of _D"(r). Substituting in Eq. ( |19[ ) we obtain 

d'^qdk I 



2d 



(27r)3 



.(G7^(q,fc)+i?j/d)|e(q,A:)r 



d^qdfc (?2|C(q, fc)|2 



(27r)3 



E 



co5(r'i(r)+V^rW) + ;.('5r(r))^ 



[l + 0(e,62,x)] 



(31) 



where C(q, /c) = X); C*; (q, fc) . 

The balance between the first |e(q, fc)p term in Eq. (|3l| ) and the 0(e) term leads to (e) ^ (p^). The O(e^) term 
depends on the distribution of e(q. A:); for thermal average it has a comparable value (section IV). 

We have identified two types of expansion parameters. The first one, x^ is related to the convergence of the I 
summation of singular vortex phases while the second one, e, is related to the response of the nonsingular phase. For 
weak Josephson coupling Aj ^ a we find x ^ (£) so that the expansion parameter is x while for Aj ^ a we find 
X ^ (e) so that the expansion parameter is (e) . 

Consider now the function (5"(r) = ^Ylii^^^)- Since q < (due to function) for p™ <C Aj the dominant 

integral range with q k. \/ p has + 77^ ~ q^. Hence (5"(r) — C"(r) — 0([p/Aj]^) and the last term in Eq. ( ^ for 
p Xj can be replaced by 



n\2 



d'vicos^l\r) ~ 1 + -Cr(r)^) = -7rln[4e](pr) 



(32) 



It is straightforward to see that the integral is convergent and therefore must be proportional to (p")^; the coefficient 
can be found after some algebra. 

The contribution to the second term in Eq. (pll) from different flux lines can be written in the form 



8 




.005 



FIG. 2. Contribution of the Josephson coupling to the energy of one displaced pancake vortex in one vortex line in units 
of 27ri?jAj (/ = Fi/(2nEjXj, see appendix C). Numerical results (dots) and the analytic form Eq. (C2) (line) are shown for 
r = 2/aV^j < 1- 



r°° rs 



^1 j ° d^qdfcg^Q(q,fc)C;(q,fc) 



n.n' l^l' d 

/— 2 
^ dkK2{vk\K - RP'|)e''=("-"')'^[l + 0(^)] 

n,n' L^f d~ 



,ik{n—n')d 

(33) 



where is the angle between p" and R" — R[5 . Since / ^ I' , the argument of the Bessel functions Ko, K2 is always 
finite and the limit in this argument can be taken; hence I ^ V terms have a harmonic expansion. 

The contribution to the second term in Eq. (|l|) from single flux lines, i.e. / ~ I' , can be computed analytically for 



/■°° n d^qdk q^\Ci{q,k)\^ 
J (27r)3 

= 2d' E P? ■ Pf / , * / J , / e^'^--'^' 



= A7rdJ2iP?f H^) + E P" ■ Pr'e-I"-"'I'''/'7I^ - n\ + Oip^X'j) . (34) 

n n=in' 

The last line is obtained by introducing x = q^JJ^ff/f and writing the integral (^) as 
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In the last two terms one can put all — > since both integrals converge. After separating the n = n' and n ^ n' 
terms and integrating over k the result Eq. (^J) is obtained. In Appendix C we consider displacement of a single 
pancake vortex in one flux line and demonstrate the agreement between a numerical exact evaluation of Eq. ( |3l| ) 
and the analytic expansion (see Fig. 2). We also show in Appendix C that for the single pancake displacement an 
expansion in p is possible also for Xj < p < a. 

We have shown in Eq. (^) that in general there are (p")^ln/9" terms in the energy expansion, confirming the 
anticipation by GK [||. Thus, strictly speaking the elastic constants are ill defined. However, Inp" is a slowly varying 
function so that replacing it by an average value In p should yield the main nonlinear correction. At finite temperatures 
p = (p^)^^^ would be a thermal average. This procedure is tested for the single pancake displacement (Appendix C) 
and is found to be in a good agreement with the exact thermal average. 



D. Effective Elastic Constants 



Effective elastic constants are obtained by replacing non-harmonic terms in p" by an average value p, which is to be 
determined self consistently e.g by a thermal average p = (p^)^/^. We introduce the effective singular phase difference 
ip"{ryff which leads to effective C'^'''^^ (r) functions 

sin^r(r)^^^ - 



[(«,"(r)2+p2)2_4(^n(j.).^n)2^2]l 

C^//(q, k) = ^ ' pqKiipq) (35) 

This approximation simplifies significantly all computations of averages because the energy can be written in the 
harmonic form 

T = T, + \j ^^{Gj\q,k)+Ej/d)\e{ci,k)-9''^ff{ci,kf 



r ' ' J (27r)3 + 4 



ttEj f^" p d^qdk 



where we define 



,da'lnHeMq,k)\' (36) 

1 +Vk 

To derive Eq. ( ^ ) the "Proper Expansion" is used with the new simJj"{rYff . The result is similar to replacing p 
by p into the coefficients of non-harmonic terms in Eq. (^l|). 

Using Eqs. (^,^ the second term in Eq. (36) can be written as 
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(fqdk |[z X q] • p{q,k)\^q^p^Kf{qp) 
X 72^" 

d ^ ' 

(Pqdk 



|[zxq].p(q,fc)|VpX'(9p) 



|p(q,fc)l' 



,1 1 

2Q^ 



(37) 



where the last integral has the analytic fit 



Qop 



1, /^pM(i±QI^ 

2^ I p-2r;2 + l 



The effective free energy of the vortex lattice can now be written in the harmonic form with effective transverse, 
C44(q, fc), and longitudinal, C44(q, A:), tilt moduli 



44(g,fc) = 044(17, fc) 



044(5, fc) = clf{q,k) + 
25$ 



2B$o 
(8^^ 
52 



In [(p-V4e)(Q2 + (i + A2,fc2)/A2)] 
1 



(SttA, 



° \n[{p^/Ae){Ql + {l + Xl,kl)/\l)\ 



(38a) 
(38b) 



It is seen that 044(9, fc) of the "naive" expansion is now corrected by replacing with p2/4e. We have assumed 
implicitly that p > ^, otherwise the "p circle" is within the vortex core area where J is reduced, i.e. the starting 
model Eq. (Q) should be modified. 



IV. THERMAL AVERAGES 



To determine the effective tilt moduli Eqs. (|3q), as well as the conditions for the expansion Eqs. 
evaluate the thermal average of the relative displacement of pancake vortices 



T 



d qdk 



sin2(f) 
q^cii + klc\^ ' q^cee + k'^c 



+ 



sin2(f ) 



we need to 



(39) 



For weak magnetic fields where Aj ^ a the Josephson contribution to C44 and 0^4 is dominant and the latter dominate 
the integrals, leading to 



4Ia5/i.(£iMM)|i + 0(4,^)1 



(40) 



where r = ^Qd/ {4tt^ X^i^) and T ^ r is assumed. The expansion parameter for Aj ^ a (Eq. ( |26b| )) together with 
Eq. ( ^ ) yields, in fact, the expansion condition T <C t. We also find by numerical integration that the displacement 
average is v? ^ p^. 

We note that p^ is nonlinear in T due to the InT factor in Eq. (ji^). Thus data on p^, e.g. by a Debye- Waller term 
in neutron scattering may probe the InT factor in weak fields, i.e. Aj ^ a. 

For strong magnetic fields where Xab ^ Xj ^ a the Josephson contribution to the last terms in Eq. (|3|) can be 
ignored, i.e. C44 



^44 



C44, leading to 



7« 1^32^[^ + ^(Af' aJ)] 



(41) 



where k = ir/d dominates the integral; here also ^ p^ . In this case the thermal average of p^ has the usual linear 
temperature dependence 

We proceed to evaluate the exp ansion parameters Eqs. (^6|) which determine the validity range of our expansion. 
For Aj » a we need x of Eq. (26a), 
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2j2 /-l/a j-njd 



X ■ 



T 



7r/d 



(Pqdk k1 

(27r)3 92(cggg2 + ct.(q^fc)fc2) 



(42) 



while for Aj ^ a we evaluate (e) directly with Eq. (^0|). Wc note that the thermal average yields (e^) = {2T/t) In A,//^ 
so that (e^) « (e). 

We find then that the eff'ective harmonic expansion is valid at temperatures below T^, 



T'^^Ta/H2Xj/Xab) 

w ra/(7rAj) 
r^«ir/ln(Aj/C) 



i/ 
i/ 
i/ 



a, Aafc ^ Aj 
a < Aj < Aafc 
Aj < a 



(43a) 
(43b) 
(43c) 



Here we defined [H jal th e decoupling temperature ~ to? ln(a/(i)/(47rA^^) for the range a, Xab ^ Aj. We note that 
the form of x Eq. ( p6a[ ) involves precisely the fluctuations that lead to the decoupling transition . It is therefore 
expected that is re lated to Td for the case when x is the relevant expansion parameter, i.e. a <C Aj. For both 
cases of Eqs. (43a, 43b) the decoupling temperature is indeed close to T** (even for a ^ Aj < Xab where the decoupfing 
transition becomes first order jj]). For a 3> Aj we expect that fluctuations of the non-singular part of the Josephson 
phase (e(q, k)) dominate so that the low temperature instability involves melting rather than decoupling. 



V. CONCLUSION 



We present in this work a proper expansion for defining elastic constants. Both deficiencies of the "naive" expansion, 
when corrected, lead to interesting physical consequences. The first difliculty is that a simple expansion at short scales, 
the "p circle" in Fig. 1, is not possible. The proper expansion shows a non- harmonic p^lnp term so that, strictly 
speaking, C44 is ill defined for displacements p > £, (for p < S, the Josephson coupling Ej should be modified by the 
reduced order parameter in the vortex core, an effect which is neglected in the Lawrence-Doniach model, Eq. (^). 

We find that effective elastic constants can be defined by replacing hip by Inp, where p is thermal average p = 
^p2^i/2^ This leads to replacing ^ of the "naive" expansion by p/\f^ in the effective C44, Eqs. (|38|). Since {p^) ^ T 
this effect can show up as a InT factor in a direct measurement of C44. Furthermore, when Xj <^ a the Josephson 
contribution dominates in the tilt moduli and {p^) ^ T/ln(r/T). This temperature dependence may be observable 
via a Debye- Waller factor in neutron scattering. 

The second deficiency of the "naive" expansion is that a Ej independent term is generated from the Josephson 
term when g — > 0. This difficulty relates to the expansion parameter x of Appendix B - the summation on flux 
lines converges only beyond a scale ~ Aj so that x ^ P^lnAj. When Ej the range where the p expansion is 
valid, X ^ 1, vanishes as 1/ In Thus at g — *■ a long range effect of many flux lines invalidates the p expansion. 
In practice one needs Aj 3> Xab for this effect to be noticeable, and the harmonic expansion is then limited to 
T <\Ta/H2Xj/Xab). 

The g — > difficulty is in fact resolved by either a self consistent harmonic approximation 0| or by a renormalization 
group method ||. In both cases the cosine function is not expanded (although ^/'"(r) is expanded as ~ 
— 2VQ!(r — R") • p") leading to a decoupling temperature Td- For T < Td Ej is renormalized to a finite value ~ 
T((^^i?j/r)^/(^~-^/^''^ which can be expanded when T <C Td/ln{T/S^'^Ej), equivalent to our expansion parameter. For 
T > Td the renormalized Ej vanishes and the p expansion is clearly invalid. 

The expansion parameter is related to Td only for Aj 3> a, while for Aj ^ a the expansion is valid for T < 
^r/ln(Aj/5). We expect that in the latter case Josephson fluxons with width Xj can form loops in between layers 
and lead to melting of the flux lattice. Thus for Xj ^ a the dominant instability is melting while for Aj ^ a it is 
decoupling; in the latter case the lattice at T > Td (held by magnetic coupling) melts at a higher temperature. 

In recent experiments on BSCCO the phase diagram has shown a number of low temperature phases (related 
to disorder |5j) while at T > AOK, where thermal fluctuations dominate, the transition to a vortex liquid phase is 
of two types: (i) At B < 500G a first order transition with no further transitions at highr temperatures, and (ii) 
at 500G < B < 900G a first order transition followed by another transition where surface barriers are reduced. For 
BSCCO Aj is estimated as 0] 800 - 2400A while at B = 500G a = 2000A. It is then consistent to consider the 
B < 500G transition as melting (Aj < a) while at i? > 500G decoupling dominates (Aj > a) with melting at a higher 
temperature. 
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A. Appendix. Tutorial Example 



Let consider two superconducting layers with Josephson coupling between them and only one pancake vortex on 
each layer. The Lawerence-Doniach free energy in the simple case of e — > has the form: 



£ = [ E [WnW]' - A72[cos(02(r) - Mr)) - 1]) 

n=l,2 



(Al) 



We now decompose the superconducting phase (?i"(r) to the nonsingular ?!>o (r) and singular part: 

Mr) =<(r)+a(r-R„) 

where a(r) = arctan(y/x) and R„ is the vortex position on the n-th layer. Define new 0{r) = [(/>2(r) — '/'i(r)] and 
V'(r) = a{r — R2) — a(r — Ri) to write the free energy in the form: 

£ = jd'r Q [V^(r)]2 - Xf [cos(0(r) + V>(r)) - 1]) + fo (A2) 

whereto part contains magnetic interaction between vortices (we do not consider it here) and i / d^r[V((/)2(r)+(/)°(r))]^ 
part which can be integrated out. 

We can perform an expansion of the cosine term with respect to the nonsingular phase 9(r) since it can be shown 
selfconsistently that |^(r)p ~ ^j^P^ ^ 1 for small relative displacements of vortices: 

p = [R2 - Ri]/2 

So we can write: 

S=lj d^r{[Ve{r)f+Xf[e{r)]'' - 2A7'[cos V'(r) - 1] + 2A7'e(r) sin V(r)) + fo + 

where 

2[p X v]^ 



sin^(r) = 



[(p2+t,2)2_4(v.p)2]i 

V = r - [Ri +R2]/2 



After shifting by the linear term in 6 term, substituting sin'^(r) by C{r) = 2[p x v]z/{p'^ + v^) we obtain: 



C(q) ,2^<Z^|C(q)P 



- Xf J d2r[cos V(r) + \{C{r)f - 1] + 0[ep^] + £0 (A3) 

Where the (?[sinV-' — C] term contributes a OiOp^) correction to the energy. 
It can be calculated analytically that: 

j d^v[cosi){T) + ^(C(r))2 - 1] = -7rln[4e]p2 (A4) 

1 f d^d g^|C(q)P 1 Xfp' , , 

2 J {2^f 1 + q^X^j - Stt^ 1 + X-/p^ + " ^ 

This shows the presence of the anharmonic term In p in the energy expansion for two superconducting layers. 
We obtain in section IIIC this anharmonicity in the more general case of a vortex lattice in a layered superconductor. 
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B. Appendix. Parameter of Expansion 



Let us introduce a 2-D lattice with I the unit ceU index and use definitions Eq. (13c) and Eq. ( [2l| ) for vp and Df-(r). 
Consider r in the I* unit cell with r = R;* + vjl , so that jv", | < a. Since (5"(r) = (5"(v") depends on r only through 
v" = r — Ri, (see Eq. (|8|)), we can write 



- V / dfc%ifi(77fe|Rz-R,.|)e 



|R/-Rp| 



x(l + 0[- 



R; — R/» 



r]) + 0(^) 



x'A = E ^e^r = I Et^TK)]' + o(4 



x(l + 0[ 



"^'wi. ft-^ IR R |J(Ri-Rp)xp,(fc)] 

-7r/d 



|R; — Re 



2 ^iknd 



R; R/* 



t]) + 0(^) 



Using the expansion 

n(i+i^r) = i + E^r + ^EE^rA' + :^EE E A"A'^r 



we obtain for the right hand side of Eq. (pfl) 



11(1 + A" W) 



(1+A") n(i+^p) 



(Bl) 



(B2) 



(B3) 



here (PviReDi — 0{p^) is used and J" means integration over the unit cell. 
For the Real part of (B2) we obtain 

I d^riRe Y[{l + Di)-1)= ( ReDi - ^ E ImDiImDi,]{l 



+ 0{[xl]\xl)) 



(B4) 



Note, that the expansion parameter depends on the configuration p". We consider averages of [x/]^ and xli which 
are diagonal in /3(q, fc), e.g. as in thermal average. We define the expansion parameter x = ([x"]^) and obtain the 
form 



Id t d'^n 

X = — / -4 dk{\p''{<i, k)\')[Maq) - Mq/Vk)f 

™^ Jbz r J-^/d 



which for Aj ^ a reduces to Eq. (26a). For Aj ^ a this is the relevant expansion parameter since from Eq. ( ^6b| ) 
X S> (e); furthermore, the other expansion parameter 



<Xr>= 271^^ In — 



(B5) 
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is seen to satisfy (x^) ^ X- 

In the case of Aj <C a the average yields 



(B6) 



which is much smaller then the other expansion parameter of Eq. (26b). Thus for Aj ^ a the relevant expansion 
parameter is (e). 



C. Appendix. One Displaced Vortex Point in One Vortex Line 

In the case of one displaced pancake vortex on the layer n = 0, in one vortex line, Z = 0, so only pj] = — = p are 
exist, we use Eq. (l27| ) in Eq. (^) to evaluate the energy numerically. The result for the energy Ti (without magnetic 
part Tv) is shown by the dots in Fig. 2. 

The numerical result can be fitted as: 

= 2nEj[1.6p' - i.oVin(-^^!ML_)] (CI) 



For p ^ A,/ we can write the energy (without magnetic part) analytically by using Eq. (|33 

= 27:Ej[L-) - p'H2pyiX'j))] (C2) 



We compare in Fig. 2 the analytic result (line) with numerical calculations (dots). 
If we put thermal average of p'^, p^ , into the non- harmonic coefRcient of p^, 

Tt" = 2nEj[1.6p' - imp' ln( y/ff ? )] (C3) 

and compute the average of p' using either effective T'l^'^ or exact JFi energy we find very similar results. Indeed, the 
approximation is good because most contribution to the average is from the region of displacements near the average 
displacement. 



[1] For a review see P.H.Kes, J. Phys. I (France) 6, 2327 (1996) 

[2] For a review see G. Blatter et ai, Rev. Mod. Phys. 66, 1125 (1995). 

[3] L. I. Glazman and A. E. Koshelev, Physica (Amsterdam) 173 C, 180 (1991). 

[4] L. L. Daemen, L. N. Bulaevskii, M. P. Maley and J. Y. Coultier, Phys. Rev. Lett. 70, 1167 (1993). 
[5] B.Horovitz and T. R. Goldin, Phys. Rev. Lett. 80, 1734 (1998). 
[6] B.Horovitz, Phys.Rev. B 47, 5947 (1992). 

[7] W. E. Lawrence and S. Doniach, in Proceedings of the Twelfth International Conference on Low Temperature Physics 

(LT-12), Kyoto, 1970, edited by E. Kanda (Keigaku, Tokyo, 1971) p. 361. 
[8] E. N. Brandt, J. Low Temp. Phys., 26, 735 (1976). 
[9] L. I. Glazman and A. E. Koshelev, Phys.Rev. B 43, 2835 (1991). 
[10] C. Kramer, Physica C 256, 236 (1996) 

[11] B. Horovitz and T. Goldin, Czech. J. Phys. 46-S3), 1833 (1996). 

[12] T. D. Fuchs, E. Zeldov, T. Tamegai, S. OOi, M. Rappaport and H. Strikman, Phys. Rev. Lett. 80, 4971 (1998). 



15 



